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ABSTRACT 

Confexf. Late-type stars interact with their close-in planets through their coronal magnetic fields. 

Aims. We introduce a theory for the interaction between the stellar and planetary fields focussing on the processes that release magnetic 
energy in the stellar coronae. 

Methods. We consider the energy dissipated by the reconnection between the stellar and planetary magnetic fields as well as that made 
available by the modulation of the magnetic helicity of the coronal field produced by the orbital motion of the planet. We estimate the 
powers released by both processes in the case of axisymmetric and non-axisymmetric, linear and non-linear force-free coronal fields 
finding that they scale as B^' 3 Bz^S^v r , where B is the mean stellar surface field, S p0 the planetary field at the poles, R p the radius of 
the planet, and v r the relative velocity between the stellar and the planetary fields. 

Results. A chromospheric hot spot or a flaring activity phased to the orbital motion of the planet are found only when the stellar field 
is axisymmetric. In the case of a non-axisymmetric field, the time modulation of the energy release is multiperiodic and can be easily 
confused with the intrinsic stellar variability. We apply our theory to the systems with some reported evidence of star-planet magnetic 
interaction finding a dissipated power at least one order of magnitude smaller than that emitted by the chromospheric hot spots. The 
phase lags between the planets and the hot spots are reproduced by our models in all the cases except for v And. 
Conclusions. The chromospheric hot spots rotating in phase with the planets cannot be explained by the energy dissipation produced 
by the interaction between stellar and planetary fields as considered by our models and require a different mechanism. 
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1. Introduction 
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Close-in planets interact with their host stars through tides and 
magnetic fields. Here we focus on the interaction between the 
coronal field of a late-type st a r and a plan et orbiting within its 
outer corona. [Shkolnik et al. (2005, 2008) presented evidence 
of a chromospheric hot spot that moves in phase with the or- 
bit of the planet rather than with the stellar rotation period in 
the cases of HD 179949 and v Andromedae. The hot spot was 
not present in all the observing seasons and its maximum vis- 
ibility was shifted with respect to the planet inferior conjunc- 
tio n. The non- steadiness of th e phenomenon was confirmed also 
by IPoppenhae ger et al. I d201 lb who did not detect any evidence 
for star-planet magnetic interaction (hereafter SPMI) when re- 
observing v And in 2009. The night-to-night variation of the 
chromospheric emission of some stars with hints of SPMI ap- 
pears to be correlated with the ratio M p sin//P or b of the pro- 
jected mass of the planet to the orbital period. This ratio can 
be assumed as a measure of the planetary magnetic moment by 
extrapolating the correlation observed in the Solar System be- 
tween the magnetic mom ent and the spin ang ular momentum 
of the magnetized bodies dShkolnik et alj |2008). Therefore, the 
proposed correlation suggests that stellar activity is somehow 
affecte d by the field of a close-by planet. However, Mille r et al.l 
d2012l) did not find evidence of activity phased with the planet 
in the case of WASP-18 whose massive (M = 10.4 Jupiter 
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masses) and close-in (orbital semimajor axis 0.02 AU) compan- 
ion is expected to induce a much larger effect than in the cases 
of HD 179949 or v And. A possible explanation could be the 
weak magnetic field of the star, as indicated by its remarkably 
low level of activity, that makes the interaction undetectable in 
spite of the supposedly strong planetary field. 

Considering a samp le of stars with planets, 
Canto Martins et al. ] d2011l) did not find any difference in 
the mean level of their chromospheric emission in comparison 
with a sample of late-type stars without detected planets, while 
iGonzalez] d201 lb claims that planet-hosting stars are, on the 
average, slightly less active than stars without planets . An 
intriguing result was recently obtained by lHartman I (120 1 ) who 
found a correlation between the chromospheric emission and 
the surface gravity of the planet in a sample of stars hosting hot 
Jupiters. 

The search for an SP MI signature in the X-r ay coronal emis- 
sion led to a controversy. Kashva p et alj ((2008) found that stars 
with close-in (orbital semi-major axis a < 0.15 AU) massive 
planets have an X-ray flux » 2 t imes higher t han stars with 
distant planets (a > 1.5 AU) and IScharf I (1201 Oh found a cor- 
relation between X-ray luminosity and the mass of the exo- 
planets, suggesting that it could be used to infer the relative 
intensity of the planetary magnetic field if the excess energy 
was released by the reconnect ion between the stellar field and 
the planetary field. However. | Poppenhaeger et al. I (l2010h and 
Popp enhaeger & Schmitt I (1201 ll) . investigating a complete sam- 
ple of planet-hosting stars within 30 pc, did not find any correla- 
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tion between stellar activity and orbital semimajor axis or planet 
mass that cannot be traced back to selection effects. 

Another possibility concerns a flaring coronal activity mod- 
ulated by the orbital motion o f the planet, as recen tly suggested 
in the case of HD 189733 bv iPillitteri etail (1201 lh . Finally, the 
possibility that SPMI manifests itself in the starspot activity ob- 
served in the phot osphere has be en also considered and is dis- 
cussed in detail bv lLanzal (1201 lh . 

Several models have been proposed t o acco unt for the pos- 
sible s ignatures of SPMI. iPreusse et al. I (120061) and [fCopp et al. 
d201 lh considered the impact on the chromosphere of Alfven 
waves excited by the orbital motion of the planet. Such waves 
can reach the surface of the star because the planet is located 
within the region where the Alfven velocity is larger than the 
speed of the stellar wind so a magnetic perturbation can move 
back to the star instead of being blown away by the wind as hap- 
pens in the case of Jupiter or Saturn in the Solar System since 
they are located in the region of the heliosphere where the solar 
wind is super-alfvenic. These models can account for the phase 
lag between the planet and a chromospheric hot spot by assum- 
ing that the spot is due to the dissipation of the waves where 
they impact onto the chromosphere, but it is very difficult to ac- 
count for the power radiated by the hot spots that is of the order 
of ~ 10 20 W (IShkolnik et al] f2005) because the wave energy is 
not focussed onto the star. Ot her models have addressed the in- 
terpretation of the phase lag (Mc lvor et al. 1120061 lLanza| |2008) 
or the mechanism of chromospheric heating (e.g., iGu & S uzuki 
l2009h . 

The problem of the energy budget of the interaction has 
been considered in nume r ical simulations, starting from the pio- 
neering work bv llp et all (120041) who, however, adopted coronal 
magnetic fields and relative orbital velocity of the planet that 
appear to be too large by about one order of magnitude than 
in the c ases of HP 179949. The sophisticated MHD simula- 
tions of Cohen et al. (2009, 2011a b) are much more realistic. 
Specifically, ICohen et al.l d201 lbh simulated the magnetic envi- 
ronment around HD 189733, a K dwarf hosting a hot Jupiter, 
approximating the large-scale variation of the radial magnetic 
field o f the star at the photosphere as derived by iFares et all 
(l2010h . They found that the power released by the star-planet 
magnetic intera ction is sufficient to a ccount for the flaring activ- 
ity detected by Pillitteri et al. (2011) or the chromospheric hot 
spots observed by Shkolnik et al. (2005) in similar systems. 

More recently. I Vidotto et al.l (1201 2l) have simulated the time 
variations of the magnetized wind of t Bootis by extrapo- 
lating the radial p h otospheric fi eld of the s tar as mapped by 
ICatala et al.l (l2007h . iDonati et alJ (l2008h . and IFares et al.1 (l2009h 
at different phases of the stellar activity cycle. The interaction 
of a strongly magnetized pre-main-sequen ce star with a close- 
in ma ssive planet has been simulated by IVidotto et al.l {2009, 
2010a). They also derived the power released by magnetic re- 
connection between the stellar and planetary fields, that could 
reach « 5 x 10 19 W for a planet orbiting at a distance of 7 - 10 
stellar radii, and estimated the expected radio emission flux and 
the effect of the strong stellar field on the migration of the planet. 

Other relevant w ork includes the po tential magnetic field 
model developed by lAdams et al.l (1201 lh to describe the mag- 
netic interaction of the components of close pre-main-sequence 
binary stars. A variant of the model has been applied to describe 
the interaction between the coronal field of a host star and the 
magn etized outflow of its close-in evaporating planet (lAdamsl 
l201lb . 

In spite of these works, a general theory of the SPMI capable 
of accounting for the required irradiated power is still lacking. A 



first step in this direction has been the work bv lLanzal (t2009) who 
found that the power dissipated by magnetic reconnection at the 
boundary surface between the planetary magnetosphere and the 
stellar coronal field is of the order of 10 17 W, i.e., insufficient by 
at least three orders of magnitude, but the interaction with the 
planetary field may trigger the dissipation of the magnetic en- 
ergy over a much larger volume by decreasing the helicity of the 
coronal field (see Sect. [2] for details). In the present paper, we 
explore this mechanism in more detail considering the magnetic 
flux tube of the stellar coronal field that interacts with the magne- 
tosphere of the planet and provide analytic formulae to estimate 
the dissipated power together with some illustrative applications 
to the systems that have shown some evidence of SPMI. Note 
that most of our considerations can be exte nded to the case of 
close binary systems with a stellar (see, e.g.. lAdams et al.ll201 ll 
iGetman et alJ l201 1 : Strassm eier et al. 11201 lh or a brown dwarf 
secondary component (e.g., iLenz et al. 1201 Oh . although here we 
shall address only the case of close-in planets. 

2. The model 

2.1. A magnetostatic coronal model 

We need a model of the corona of the host star to test the validity 
of our assumptions on the energetics of the star-planet interac- 
tion. Assuming that the coronal plasma is in hydrostatic equilib- 
rium under the action of the pressure gradient, the gravity, and 
the Lorentz force, we have: 



Vp + pV$ + J x B - 0, 



(D 



where p is the plasma pressure, p its density, O the total poten- 
tial (gravitational plus centrifugal), B — Bs the magnetic field, 
with s the unit vector in the direction of the magnetic field, and 
J the current density. The centrifugal potential can be neglected 
in our case, provided that the star has a rotation period of at 
least ~ 10 days and a radius and mass comparable with those of 
the Sun. Therefore, O =s GM/r, where G is the gravitation con- 
stant, M the mass of the star, and r the distance from its centre. 
The thermal conduction is high along magnetic field lin es and 
strong ly inhibited in the orthogonal direction (cf., e.g., Priest 
1984), thus we assume that the temperature T is constant along 
magnetic field lines, i.e., dT/ds = 0. The dot product of Eq. ([1) 
by the unit vector s gives: 



dp dO f . 
os or 



0. 



Since dr — (r ■ s)ds, this becomes: 

dp (9$ n 
or or 



(2) 



(3) 



that can be integrated along a given magnetic field line by con- 
sidering the ideal gas law to eliminate p - pp/RT, where p is 
the mean molecular weight of the coronal plasma, R the gas con- 
stant, and T the temperature along the given field line. Thus, the 
variation of the pressure along a field line is: 



p(r) = p exp 



(4) 



where R is the radius at base of the corona, that we assume to 
coincide with the radius of the star, and Hq = (RTR 2 )/{GpM) 
the pressure scale height at the base of the corona. For the Sun, 
p = 0.6 and H = 5.1 X 10 7 r m, with T in MK. The variation 
of the density can be derived from the ideal gas law and has the 
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same radial dependence of the pressure along a given field line. 
Note, however, that T, po, and the base density po vary in general 
from one field line to the other. We shall apply this simple model 
in Sect. 13. 1 I to justify the assumptions of our energy dissipation 
model. 



2.2. Power dissipated by magnetic reconnection 

The orbital motion of the planet inside the stellar corona pro- 
duces a continuous reconnection betw een the coron al magnetic 
field lines and the planetary field lines. lLanzal (120091) studied this 
phenomenon and estimated the power dissipated in the stellar 
corona. We briefly recall the main assumptions and the results 
of that investigation that are useful for the present study. We as- 
sume that the planet is on a circular orbit located on the equa- 
torial plane of the star. The surface of the planetary magneto- 
sphere, where its field lines interact with those of the coronal 
field, is assumed to be a sphere of radius R m . As a matter of 
fact, the magnetospheric boundary can be elongat ed in the di- 
rection of the orbital motion of the planet (cf., e.g., ICohen et al.l 
1201 However, we specialize our theory to the case when the 
orbital velocity of the planet is much smaller than the Alfven 
velocity in the stellar corona. Therefore, the magnetic field con- 
figuration can be regarded as magnetostatic and the boundary of 
the magnetosphere is nearly spherical because it is defined by 
the balance between the magnetic pressure of the coronal field 
and the pressure of the planetary field assumed to be a dipole. 
Assuming that the magnetic pressure of the coronal field B(r m ) 
at the boundary of the magnetosphere is in equilibrium with the 
pressure of the planetary field B p (r m ), we have: 



B 2 p (r m ) = B 2 (r m ), 



(5) 



where r m is the position vector of a generic point on the bound- 
ary of the magnetosphere. The planetary field can be assumed to 
be that of a dipole, so that its variation with the distance A from 
the centre of the planet is: 



-3 



(6) 



where B p q is the field at the poles of the planet and R p its radius. 
At the boundary of the magnetosphere, we have A = R m and, 
considering Eq. Q, we find: 



Rm = Rr 



B(r p ) 



»p0 



-1/3 



(7) 



where B(r p ) is the intensity of the ambient coronal field at the 
position r p of the planet. Since the radius of the magnetosphere 
is generally small in comparison with the lengthscale of vari- 
ation of B in the outer stellar corona, we neglect the variation 
of B across the magnetosphere and consider its intensity at the 
location of the planet r p . 

The power released by the magnetic reconnection at the 
boundary of the planetary magnetosphere can be estimated as: 



Trec-B (r p )R m v rd 

ft 



y W c-R 2 p B AI \r p )B 2 ^v KX , 



n 

.-7 
A* 



(8) 



where ju is the magnetic permeability of the plasma, < y KC < 1 
an efficiency factor that depends o n the angle b etween the inter- 
acting magnetic field lines (e. g., IPriest Il2003l) . v re i the relative 



velocity between the reconnecting magnetic field lines, i.e., the 
planetary and the stellar fields, and we have assumed that an 
effective surface nR^ is available for the interaction of the re- 
connecting field lines. 

The powe r dissipated i n the reconnection process has been 
estimated by lLanzal d2009t) to be of the order of 10 17 W in the 
case of a stellar dipolar potential field, i.e., insufficient by ~ 3 
orders of magnitude. In Sect. [3] we shall see that considering 
force-free magnetic fields we can increase that power by about 
two orders of magnitude in the most extreme cases. This is still 
insufficient to account for the power irradiated by the chromo- 
spheric hot spots. Therefore, we shall consider another mecha- 
nism that can release greater powers by extracting energy from 
a larger coronal volume, not only from that taking part in the re- 
connection. This mechanism is connected with the role of mag- 
netic helicity in stellar coronae that we shall briefly describe in 
the next Section. 



2.3. The role of magnetic helicity in stellar coronae 

We assume that the Lorentz force is dominating over all the 
other forces in the region of the stellar corona where a close- 
in planet is located. This is valid provided that the ratio between 
the plasma pressure and the magnetic pressure B = 2/jp/B 2 is 
much smaller than the unity in the considered coronal domain, 
as we shall show in Sect. 13.11 With this assumption, the force- 
free approximation can be applied to describe the coronal field, 
i. e., 



VxB = aB, 



(9) 



where the force-free parameter a is constant along each field 
line, as immediately follows from the curl of the defining equa- 
tion (0. In general, a will vary from one field line to the other 
(non-linear force-free field). When it is constant for all the field 
lines, the field is said to be a linear force-free field. We assume 
that the orbital motion of the planet and the associated plasma 
flow do not significantly perturb the stellar field configuration 
which is determined by the boundary conditions set at the base 
of the corona on the stellar surface. If those boundary condi- 
tions evolve on a timescale much longer than the Alfven travel 
time across the corona, the field can be assumed to be at each 
instant in a magnetostatic configuration as described by Eq. (0. 
Of course, the force-free approximation is not valid at the stel- 
lar photosphere, where the pressure of the plasma is comparable 
or greater than the magnetic pressure, but we can assume that 
our low-beta approximation is valid starting from the base of the 
corona that we assume for simplicity to coincide with the stel- 
lar surface because the photospheric and chromospheric pres- 
sure scale heights are much smaller than the radius of the star. 
Note also that the flux systems of the stellar corona and of the 
planetary magnetosphere are topologically separate because the 
planetary magnetic field is potential close to the surface of the 
planet. Therefore, a stationary magnetic field line that intercon- 
nects the planetary field with the stellar coronal field would have 
a zero value of a that is in general incompatible with the pres- 
ence of electric currents flowing through the stellar corona. In 
other words, the field lines of the planetary field must be con- 
fined within the planetary magnetosphere and interact with the 
stellar field lines only on the boundary of the magnetosphere 
where a time-dependent reconnection occurs. A large-scale flux 
system interconnecting the stellar and planetary fields in a low- 
beta regime is possible only when the stella r corona l field i s as- 
sumed to be poten tial as in the models of lAdamsl (1201 ll) and 
lAdams et al.1d201lh . 
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If the magnetic field is confined to some closed volume V, 
i.e., its field lines do not cross the boundary S of V, it is possible 
to define a conserved quantity in ideal MHD called magnetic 
helicity H: 

H= \ A-BdV, (10) 

where B = V x A is the magnetic field and A its vector potential. 
H is a topological measure of the twisting of th e magnetic field 
lines and of their degree of cross linkage (cf. iBerger & F ield 
Il984t Demoulin et al. 2006, for more details). The conserva- 
tion of the helicity means that, if the magnetic Reynolds number 
is very large, so that the diffusion of the field is negligible in 
comparison with the induction effects, the value of H does not 
change vs. time whatever the plasma motions inside the volume 
V. Moreover, the minimum energy state allowed in the case of a 
finite helicity is a linear force-free field dWoltierlll958l) . instead 
of a potential field which would represent the absolute mini- 
mum that c an be reached only when the helicity is zero (cf., e.g., 
iPriest 111984 Ch. 3). 

In a non-ideal plasma, the finite resistivity produces a decay 
of the field with a conversion of the magnetic energy into ther- 
mal and kinetic energies. Experiments with laboratory plasma 
and theoretical considerations have shown that the variation of 
the magnetic helicity during the relaxation to the minimum en- 
ergy state is extremely sl ow in compari son with the magnetic 
energy dissipation (e.g., IBerger ill 9841) ; in practice, the relax- 
ation proceeds under the constraint determined by the helicity 
conservation and the final relaxed state is a linear f orce-free field 
with t he same amount of helicity as the initial field dTavlorll974L 
[19861) . 

The general results recalled above can be applied to a stel- 
lar corona if we account for the differences between a laboratory 
plasma and a coronal plasma. Firstly, the coronal field lines are 
not confined by a magnetic surface, i.e., a surface over which the 
normal component of the field vanishes. They cross the surface 
of the star so that the magnetic helicity is no longer a gauge in- 
variant quantity, i.e., it varies by adding the gradient of a generic 
scalar function^ to the vector A leading to the ga uge transforma- 
tion A — > A + Vy. To ove rcome this difficulty, Berger & Field 
(119841) and IBerger I (1 1985b defined a relative magnetic helicity 
that is independent of the gauge assumed for the vector potential 
and is conserved in ideal MHD. In a finite domain, the mini- 
mum energy field having a given rel ative helicity i s again a lin- 
ear force-free field (cf. Sect. Ill of IBerger |["l985l) . The second 
difference between a laboratory plasma and a coronal plasma is 
that the coronal field can extend to the infinity. In this case the 
minimum energy state is actually reached by driving the twist 
of the field lines to the infinity which dilutes the magnetic he- 
licity density until the final state is virtually indistinguishable 
from the potential field corresponding to the boundary condi- 
tions applied at th e coronal base (cf., e.g., IVekstein et al.lll993t 
IZhang et~a"D l2006j and references therein). However, a stellar 
corona consists also of closed magnetic structures that are con- 
fined by overlying open fields, as in the case of a magnetic arcade 
with an helmet streamer on its top. In those closed structures, 
the accumulation of magnetic energy due to the shearing mo- 
tions of the photospheric footpoints of the field lines is accom- 
panied by a relaxation produced by the dissipation of magnetic 
energy localized in thin current sheets. The minimum energy 
state of each confined structure is predicted to be a linear force- 
free field satisfying the con straint of relative helicity conserva- 
tion (see iDixon etaill 19891 for details). However, if the amount 
of helicity accumulated in a confined structure exceeds a certain 



threshold, that depends on the boundary conditions, the mag- 
netic configuration may become unstable and erupts producing 
a coronal mass ejection (hereafter CME) that takes away most 
of the accumulated helicity al lowing the remaining field to relax 
to a quasi-potential state (cf . IZhang et ai1l2006l : IZhan g & Fiver 
l2008tlMiller eTaIll2009T) . 

The processes that dissipate the excess magnetic energy and 
lead to the linear force-free minimum energy state are outside of 
the scope of the theory described above. In general, they involve 
the generation of hydromagnetic instabilities and/or the forma- 
tion of localized currents under the action of photos pheric mo- 
tions or the emergence of new magnetic flux (e.g., Pries t et ail 
120051 : iBrowning et al. 1 120081) . Since the timescale for attain- 
ing that minimum energy state is of pa rticular relevance to our 
model, we briefly refer to the work of IBrowning et al. I ([2008 ) 
who suggest that the non-linear developments of ideal magneto- 
hydrodynamic instabilities are the best candidates for the dissi- 
pation of the excess energy in confined structures. The reason for 
choosing ideal rather than resistive instabilities is that only the 
former have sufficiently short timescales in the highly conduct- 
ing coronal plasma to be relevant to flares or coronal heating. 
Their 3D numerical simulations of the non-linear evolution of 
a kink instability in a magnetic flux tube show that the field is 
driven toward the minimum energy state while its relative he- 
licity is conserved with a characteristic time scale comparable 
with the Alfven transit time from one base of the flux tube to the 
other. This compares well with the observations of solar flares 
showing that most of the energy released in the impulsive phase 
is dissipated over timescales of the order of 10 - 100 s in flux 
tubes having lengths ranging from 10 7 to 10 8 m, i.e., their relax- 
ation occurs on timescales comparable with the Alfven transit 
time along the field lines since the typical Alfven velocity in the 
solar corona is of the order of 10 6 m s . 



2.4. Assumptions on the characteristic timescales 

We shall consider the variation of the relative helicity produced 
by the orbital motion of the planet across the coronal field of 
its host star and estimate the energy made available for dissi- 
pation by this process. The application of the above theory is 
considerably simplified when the timescale tu for the helicity 
variation is considerably longer than the Alfven travel time t\ 
along the magnetic field lines from the star to the planet. On 
the other hand, the timescale ?r for the relaxation to the mini- 
mum energy state under the constraint of helicity conservation 
can be assumed to be comparable with t& (cf. the final paragraph 
of Sect. 12.3b . Therefore, we can assume that the coronal field is 
always close to the minimum energy state determined by its in- 
stantaneous relative helicity. These hypotheses allow us to study 
the evolution of the coronal field as a sequence of magnetostatic 
configurations to which a simple model for the energy dissipa- 
tion can be applied (see Sects. 12. 5 land l2.6l l. 



2.5. Connection between helicity and magnetic energy 

The relative helicity of a coronal magnetic field configu- 
ration can be comput ed following the method described in 
Berger & Field I (11984 . Specifically, we consider a magnetic 
flux tube connecting the surface of the star with the magneto- 
sphere of the planet (cf. Fig. [TJ. Its bases are F\ on the surface 
and F2 on the magnetospheric boundary. Using the formulation 
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of lDemoulin etalJ (2006), its relative magnetic helicity can be 
written as: 



(11) 



H R = f (A + A p ) • (B - B p ) dV, 

Jv, 

where Vf is the volume of the magnetic flux tube, B = VxA the 
magnetic field, A its vector potential, B p = V x A p the potential 
magnetic field having the same normal component of B on the 
boundary of Vf, and A p its vector potential (see Appendix IA1. 
We adopt the so-called Coulomb gauge, i.e., V ■ A = and V ■ 
A p = 0. Applying vector calculus identities and Gauss' theorem, 
Eq. ( fTTT i can be recast as: 

H R = I A ■ B dV + (h (Ax A p - ifrAp) ■ h dS, (12) 
JV, Js, 

where Sf is the closed surface bounding the volume Vf, h the 
unit outward normal to Sf, and iff the scalar potential of B p , i.e., 
B p =Vi//. 

The connection between the magnetic energy and the relative 
helicity can be derived by considering that V ■ (A x B) — B ■ (V X 
A)-A-(VxB) = B 2 -a(A-B), where we have applied the force- 
free condition V x B = aB to the field of the flux tube. Since a 
varies only across the field lines, the variation of a across the flux 
tube section is limited, i.e., a m \ n < a < a m&x , with the minimum 
and maximum values close to each other if the cross section is 
sufficiently small. The magnetic energy E of the flux tube is: 



E = \ T- = T- f a( - A - B )dV + 
>V, 2 V 2 V Jv, 



- f 

2 M Jv, 



V -(AxB)dV 



= — f a(A ■ B) dV + — £ (A x B) ■ MS ,(13) 

2^ JV, 2 V Js, 

where the second equality follows by the application of Gauss' 
theorem. Applying the mean-value theorem to the first integral, 
we find: 



/ 

Jv, 



—H R + - 

2fi 2fi 



ABdV + ^- £ (AxB) -MS 

2 V Js, 



[{a)(ifrA p - Ax A p ) + AxB] - MS, (14) 



where the mean value a m [ n < {a) < a m - dx and we have made 
use of Eq. (fT2l to transform the fir st integral on th e r.h.s. This 
expression generalizes Eq. (17) of iBergerl (119851) to the case 
of a non-linear force-free field. Since the magnetic field inside 
the flux tube occupies a finite volume and is confined by the 
closed surface Sf, its minimum energy state is the linear force- 
free field satisfying the boundary conditions and the constraint 
set by the conservation of the relative helicity. If the dissipation 
of the excess magnetic energy is fast in comparison with the he- 
licity variation, we can assume that the field is always close to 
such a minimum energy state. A variation of the relative helicity 
or of the boundary conditions that changes the surface integral 
in the r.h.s. of Eq. (fT~4T > then produces a variation of the energy of 
the field. The contribution of the surface integral to the energy 
can be written in terms of the two quantities: 



= 



(f) (i/'Ap - Ax Ap) ■ hdS, and 

Js, 



2 = — £ (AxB) - MS, 

2 HJs, 



(15) 
(16) 



We focus on the effects produced by the orbital motion of the 
planet that changes the boundary conditions on the base F%, The 



variation of the surface terms is confined to F2 because we as- 
sume that the stellar coronal field is not perturbed outside the 
planetary magnetosphere by the motion of the planet. This gives: 



dE _ (a) (dH R d<3>\ dL 
dt 2yU \ dt dt J dt' 



(17) 



where (a) is assumed to be independent of the time because a is 
constant along each magnetic field line so that its mean value is 
set by the boundary conditions at F\ on the stellar surface that is 
not perturbed by the planet. For the sake of simplicity, we com- 
pute the variati ons of Q and X by assuming that the surface F2 is 
fixed (see, e.g., ISmirnov 1[T%1 for the neglected terms). As A p 
and \jj are completely determined by the normal component of 
the magnetic field B n that does not change versus the time (see 
Sect. 12.6b . their time derivatives vanish. Therefore we have: 



(18) 



d® r 

— = - [(v x B)x A ? ] hdS = 

dt J/7, 

= JT [(A p -B)v n -(A p -v)B n ] dS, 

where we have made use of ^ = v x B in the Coulomb gauge 
and v n and B a are the normal components of the velocity and the 
magnetic field on F^. Similarly, the variation of the surface term 
2 is: 



dl, 1 d r 
dt 2fi dt J F 

2 M Jf 2 ^ 



(AxB)- MS = 
xB)xB + Ax[Vx(vxB)])-MS, (19) 



where we have made use of the induction equation of ideal MHD 
and of ^ = v x B. The contribution of dZ/dt to the energy 
variation is therefore of the order of vB 2 F2l2ji. Considering a 
maximum cross section of the flux tube equal to the cross section 
of the planetary magnetosphere, i.e., F2 = ttR^, the contribution 
to the power is: 



cTL 
dt 



■2 d2 , 



2^ 



-B A R 



(20) 



This is comparable with the energy released by the reconnection 
between the stellar and planetary magnetic fields on the bound- 
ary of the planetary magnetosphere (cf. Eq. |8), as discussed in 
Sect.O 

Considering this estimate of the magnitude of d~Lfdt, we see 
that when 2(a)(A p ■ B) » B 2 the terms containing the derivatives 
of the magnetic helicity and of dominate over the derivative of 
X in Eq. (TTTl i and we can write: 



dE^ia^ldHR d® 
dt 2a \ dt dt 



(21) 



Note that Eq. (I2TI) provides a lower limit for the dissipated power 
because it was obtained in the hypothesis that the field is always 
close to the minimum energy state corresponding to its relative 
helicity. If the field has a larger excess energy, i.e., it is in a non- 
linear force-free state far from the minimum energy state, then a 
greater energy can be released. 

2.6. Variation of the magnetic helicity and energy 

The relative magnetic helicity varies according to the follow- 
ing formula derived by applying the induction equation of ideal 
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Fig. 1. A sketch of the star-planet system with indication of V a , the volume interior to the star, and V m , the volume occupied by the 
planetary magnetosphere assumed to be spherical. Two flux tubes are plotted and in the upper one Vf the cross sections F[ and Fi 
where the tube intersects the stellar surface and the magnetosphere are labelled, respectively. The blue arrows indicate the relative 
orbital velocity v and its normal component v n at the intersections of the two flux tubes with the magnetosphere. The dashed line 
inside V m separates the region I where v n < 0, i.e., is oppositely directed with respect to the outward normal to the base of the flux 
tube, from the region II where v„ > 0. 



MHD to the definition of relative helicity (see Berger & Field 
1984 lVeksteinetaD[T99l : 



dH R 
dt 



-2^[(A p -B)v n -(A p -v)B n ] dS, 



(22) 



where S is the boundary of the domain occupied by the field and 
v n and B n the normal components of the plasma flow v and of 
the magnetic field B on S , and A p the vector potential of the 
potential magnetic field with the same normal component as B 
on S ; the unit normal is taken positive in the direction outward 
the volume occupied by the field. 

The emergence of the magnetic flux from the convection 
zone of the star and the velocity fields at the base of the corona 
produce a continuous variation of the relative helicity that can be 
associated with the energy dissipation responsible for the coro- 
nal heating and the storage of excess energy to be delivered in 
CMEs. These processes are independent of the presence of the 
planet and contribute to the unperturbed energy budget of the 
stellar corona. Here we focus on the processes related to the 
planet. In the framework of our hypotheses, they are related to 
the magnetic and the velocity fields on the interface S m between 
the planetary magnetosphere and the coronal field. The surface 



S m is a closed surface and occupies a limited volume of the 
corona. For the sake of simplicity, we have assumed it to have 
a spherical shape. The magnetic field is tangent to S m because 
it separates the coronal field from the planetary field and the 
pressure of the latter forbids the penetration of the former into 
the planetary magnetosphere. Since the radius of S m is small in 
comparison with the lengthscale across which the coronal field 
varies, we can assume that the local field at a generic point on S m 
is B m - B + B', where B is the unperturbed field of the corona 
(assumed uniform) and B' a local perturbation that makes the 
local field tangent to the magnetosphere, i.e., B m ■ h s = 0, where 
n s is the normal to S m . Since the magnetic energy is steadily dis- 
sipated by reconnection on the boundary of the magnetosphere, 
the energy of the field B m is minimized. The minimum energy of 
the field B m compatible with the prescribed boundary condition 
is obtained when B' is a potential field, i.e., B' = V£, where the 
potential £ satisfies Laplace equation V 2 £ = and B' decreases 
to zero at the infinity since the effects of the magnetosphere are 
localized. The solutions of Laplace equation that vanish at the 
infinity are proportional to 1 jr and its derivatives with respect to 
the spatial coordinates, where r is the distance from the centre of 
the planet. Given the complete symmetry of a spherical magne- 
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tosphere, only the constant vector B can appear in the solution 
and, considering the linearity of both Laplace equation and the 
boundary conditions, f must involve B linearly. The only scalar 
that can be formed from B and the derivatives of 1 jr is the scalar 
product B ■ V(l/r) (cf. the mathematically analogous case of the 
potential flow of an ideal incompressib le fluid around a spherica l 
body moving through the fluid, e.g., lLandau & LifshitzJll959h . 
Therefore, we seek f in the form: 



because B„ 



1 



? = CV - =-(C-r)-, 



where C is a constant vector that is chosen to satisfy the bound- 
ary condition on the sphere S m , i.e., B' r — -B ■ r. Since 
B r = d£/dr, we immediately find C = — (7?^/2)B yielding: 



2 r L 
The perturbation of the magnetic field is: 



B' 



V ^ 2 U™ 



[B - 3(B ■ r)r] 



(24) 



(25) 



The variation of the magnetic energy of the coronal field pro- 
duced by the perturbation B' is (see Appendix IBT>: 



AE 



~V m B 2 
An 



-—R p B p qB, 



(26) 



where V m = (47r/3)/?„ is the volume of the magnetosphere, R p 
the radius of the planet, B p q the field intensity at the pole of the 
planet, and B the intensity of the coronal field at the location of 
the planet. Considering a magnetic loop that is crossed by the 
planetary magnetosphere, its energy will change by AE over a 
timescale r cross = 2J? m /v re i, where v re i is the relative velocity be- 
tween the planet and the field. Therefore, the power dissipated 



is: P 



\AE\/t c 



{l/2)(7i/fj.)Rf n B Vrei, i.e., it is the same 



as the power dissipated by magnetic reconnection at the mag- 
netospheric boundary for y lec = 1/2 (cf. Eq. [8]). As we saw in 
Sect. 12.21 it is insufficient for our purposes, thus we proceed with 
evaluating the energy released by the magnetic helicity variation. 

The potential magnetic field B pm satisfying the same bound- 
ary condition of B m on S m can be found with a similar argument. 
We consider its vector potential A pm = A p + A p , where A p is the 
vector potential of the unperturbed potential field B p that de- 
pends on the boundary conditions at the stellar surface, and A' p 
is its local perturbation. The local perturbation of the potential 
field B' p = V x A' p can be derived by the same argument applied 
above to find B' yielding B' p = (l/2)(r/,ff m r 3 [B p - (B p ■ r)r}. 
The corresponding vector potential in the Coulomb gauge is: 



A' 



-{- 
2\R m 



(r x B p ). 



(27) 



As in Sect. 12.51 we consider a magnetic flux tube connecting the 
surface of the star with the magnetospheric boundary S m and 
indicate by F2 its base on S m (cf. Fig.[T|). The rate of change of 
the magnetic energy of the flux tube produced by the variation of 
the relative helicity and of the surface term can be computed 
by Eqs. ([21), 423, and ([T8j and is: 



= § F [ (Ap + A ' p) ' (B + V " dS ' 



(28) 



-B m h s — on the boundary of the magneto- 



sphere. Since B and A p can be considered uniform on S m , on 
account of Eqs. ( f25l ) and ( l27l i. we obtain: 



dE 
dt 



~ \[{rxB p )-B]v n dS. 



(29) 



(23) The second integral on the r.h.s. can be neglected in comparison 



to the first one because \r x B n 



R m B p is much smaller than 



|A p | that is of the order of aB p , where a is the semimajor axis of 
the planetary orbit. 

To compute the first integral, we adopt a Cartesian reference 
frame with the z axis along the relative orbital velocity between 
the planet and the stellar field v and the x axis chosen so that 
the vector B lies in the xz plane. Moreover, we choose the base 
F2 of the connecting flux tube to coincide with the hemisphere 
of S m where v n has a constant, say, negative, sign (cf. Fig. [TJ. 
Performing the integration over that hemisphere, we find: 



dE 
dt 



3n(a) 



1 



— = -^R z m v re] [(A p ■ B) + -(2B z A pz - B x A px )], 



(30) 



where v re i is the relative velocity between the planet and the stel- 
lar coronal field. When the stellar field is axisymmetric, the orbit 
of the planet is circular and lying in the equatorial plane of the 
star, the vector potential A p = A p z (cf. AppendixlAl) and Eq. (f30b 
simplifies to: 



dE 
dt 



9n(a) , 

-A^v rel (A p . B). 



(31) 



2.7. Application to the star-planet interaction 



Two kinds of interaction occur on the boundary of the planetary 
magnetosphere S m from the point of view of the helicity varia- 
tion. Specifically, a coronal flux tube whose field lines touch the 
magnetosphere in the domain where v n < 0, i.e., on the left of 
the dashed line inside V m in Fig. [1] will experience an increase 
of its relative helicity which reduces the amount of magnetic free 
energy, thus opposing dissipation. This helicity increase can trig- 
ger a CME with an associated flare if the helicity previously ac- 
cumulated is close to the thresh old for the loss of equilibrium of 
the field configuration (cf., e.g., IZhang et al.l2006l) . On the other 
hand, a flux tube whose field lines touch the magnetosphere in 
the region where v n > (i.e., on the right of the dashed line in 
Fig. [TJ will experience a decrease of the relative magnetic he- 
licity leading to an increase of the magnetic free energy and an 
enhancement of its dissipation. In the case of a coronal loop as 
sketched in Fig. \T\ the helicity increase and decrease in the two 
legs compensate for each other because the magnetic field re- 
covers the initial unperturbed configuration after the passage of 
the planet through the top of the loop. The only net effect is a 
modulation of the magnetic energy dissipation that is initially 
reduced and then increased by the passage of the planetary mag- 
netosphere across the coronal field lines. The total energy budget 
of the corona is not affected by this process, thus the X-ray lumi- 
nosity of the star is not affected when averaged along one orbital 
period of the planet. Note, however, that the increase of helic- 
ity produced in the first part of the modulation, i.e., in the flux 
tubes for which v„ < 0, may trigger a CME or a flare when a 
previous accumulation of helicity has brought the magnetic con- 
figuration close to the threshold for instability. This mechanism 
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may in princi ple explain the mod ulation of the flaring activity 
suggested bv lPillitteri et all (1201 ll) in the case of HD 189733. 

These predictions are valid only if the timescale for the 
energy dissipation ?r is remarkably shorter than the timescale 
for the helicity variation induced by the planetary motion, i.e., 
tR <s tu, otherwise the effects of the helicity increase and de- 
crease would be averaged to zero along the timescale ?r and no 
significant energy dissipation could be observed. Moreover, if 
tR > tu, the phase lag between the planet and the chromospheric 
hot spot should continuously vary because of the statistical char- 
acter of the energy dissipation that implies a range of delays 
between the perturbation induced by the planet and the energy 
release. 

In our treatment, we have considered only the velocity field 
arising from the orbital motion of the planet across the magnetic 
field of the stellar corona. However, other velocity fields may 
potentially be relevant for the variation of the magnetic helic- 
ity, e. g., those associated with the reconnection of the magnetic 
field lines or the evaporation of the planetary atmosphere. The 
reconnected field lines must leave the region where reconnec- 
tion has occurred and this produces a flow that is nearly orthog- 
onal to the incoming flow, i.e., the flow carrying the opposite 
field lines into the reconnection region. Its speed is of the order 
of the local Alfven speed, but in our case it is nearly tangent to 
the magnetospheric boundary with a small component outward 
from it because the magnetic pressure of the planetary field halts 
the penetration of the reconnected field lines into the magneto- 
sphere. Therefore, this velocity field will lead to v n < produc- 
ing an increase of helicity. The same is true for the evaporation 
flow of the planetary atmosphere that has been observed in some 
systems, notably HD 209458 and HD 189733, and is thought 
to be induced by the strong irradiation by the close host star 
(e. g., Lecavelier des Etangs et al.ll2004t lEhrenreich et al.1 12008: 



iLinskv et al.ll2010l) . The speed of the evaporation flow at the 
limit of the magnetospheric boundary is still uncertain, but a 
reasonable estimate is a few tens of km s _1 , i.e., significa ntly 
smaller than the orbital velocity of the planet (Adams 20 lT]). In 
conclusion, the main contribution to the modulation of the mag- 
netic helicity comes from the orbital motion. 

2.8. Magnetic field configurations considered for the 
computation of the energy dissipation rate 

To compute the amplitude of the modulation of the energy dis- 
sipation rate by means of Eqs. (l30l l and (IBTb . we now assume 
a spherical polar reference frame with coordinates (r,0,cp), the 
origin at the centre of the star and the polar axis along the stellar 
rotation axis. Our reference frame rotates with the stellar angular 
velocity Q with respect to a distant observer. 

To obtain quantitative predictions from our model, we need 
to specify the configuration of the stellar magnetic field. Linear 
force-free fields obey an Helmoltz equation and can be expressed 
analytically both in the axisymmetric and non-axisymmetric 
cases allowing us to study the role of axisymmetry on the SPMI. 
Their main limitation is that they consist of a set of infinite dis- 
joint subdomains each of which is confined between two con- 
centric spherical surfaces. In other words, we must restrict our 
consideration to a field between the su rface of the sta r and some 
magnetic surface at a radius as in lLanzal <|2009). This field 
cannot extend to the infinity which leads to a severe restriction 
imposed on the topology of the coronal field. Nevertheless, in 
view of their mathematical simplicity and the possibility of treat- 
ing non-axisymmetric configurations, we shall consider in detail 
linear force-free fields. On the other hand, non-linear force-free 



fields pose formidable mathematical problems. Therefore, we 
shall limit ourselve s to the simple configurations introduced by 
lLow & Lou 1 (11990) and then consider an extension o f our results 
to the more general configurations introduced by IFlver et aTTI 
(2004). 



2.9. Linear force-free fields 

First, we consider a linear forc e-free field as introduced by 
IChand rasekhar & Kendaii1 (ll957l) . The dependence of the vector 



potential A p on the distance r from the star is given by Eqs. (IA.4t 
and iA.li in Appendix lAl It varies as ir/R)~" where n > 1 is 
the order of the field multipole and R is the radius of the star. 
Since r/R ~ 7 - 10 in the case of close-in planets, the leading 
term is that corresponding to the dipole, i.e., n = 1, and we can 
neglect the terms of order n > 1 because their relative contribu- 
tions to the energy dissipation rate are of the order of ~ (r/Ry 2n . 
With these assumptions and approximations, we can estimate 
the order of magnitude of the dissipated power by considering 
for simplicity only the term A p • B at the location of the planet 
r p = (r,9,<p) in Eqs. (l30l and (T3TT >; we find (cf. Appendixes lAl 
and0: 



-2 



T — ( — ) 8(l) \ B l sin2 f - BiBo sin 26>cos(0 - f) • 
4 a \ q j y 

0]}. 



£) + cos 2 0cos 2 ( 



(32) 



+ B 2 [sin : 

where Bo is the intensity of the axisymmetric component of the 
magnetic field (i.e., the mode with m - 0) at the poles of the star, 
B\ the intensity of the radial component of the non-axisymmetric 
field (i.e., the mode with m = 1) on the equatorial plane (0 = |) 
at the longitude a the force-free parameter, qo = \a\R, q = \a\r, 
and the func tion g(q) has been defined in Eq. (2) in Sect. 3 of 
lLanzal d2009l) . For completeness, we provide the expressions of 
the magnetic field components in AppendixICltogether with that 
of g(q). 

Assuming that the planetary orbit is on the equatorial plane 
(6* = 7r/2) and is circular with a radius a, the expression of (A p -B) 
simplifies to: 

(A p ■ B) rp = l -aR 2 g(g) [Bl + B\ sin 2 (0 - £)] , (33) 

where q = \a\a. 

The squared intensity of the stellar magnetic field B at the 
distance of the planet is given by (cf. AppendixICb: 



B\r v ) = - 



1 


(4) 


4 






\ 4 




? ' 



[lg'(q)] 2 + lg(q)] 2 } [B 2 + B\ sin 2 (0 - £)] 



l2 D 2, 



(34) 



The magnetospheric radius R m follows by substituting Eq. (1341 
into Eq. Q. We immediately see that R m is not constant when 
Bi + 0, but is modulated by the terms cos 2 (0-^) and sin 2 (^>-£), 
i.e., it has a time frequency of 2(w - Q), where to is the orbital 
frequency of the planet. 

Eqs. d33l and (0 provide the basic ingredients to compute 
the energy dissipation rate from Eqs. (l30l and (13TV It is useful 
to consider separately the cases of an axisymmetric coronal field 
and that of a non-axisymmetric field to point out the difference 
in the frequency of the modulation of the magnetic dissipation 
rate. 
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2.9.1. Axisymmetric linear force-free field 

When |Bq| » \Bi\, the axisymmetric field dominates and the am- 
plitude of the modulation of the energy dissipation rate as given 
by Eq. OTb becomes: 

Phel = ^4-i-R^BlBl Q qJq-ig(q){[g'(qyf + [giq)] 2 }^ V rel .(35) 

This expression has the same dependence on R p , Z?o, B p o, and 
v re i of the energy dissipation rate produced by the reconnection 
between the coronal and the planetary fields (cf. Sect, \2.2l . The 
ratio between the two powers is: 



fhel 9 ^ 



g(q) 



8y.ec \ q ) [g'(q)] 2 + [g(q)] 



2' 



(36) 



where q = \a\a is the non-dimensional orbital radius of the 
planet. In the distance range where close-in planets are usually 
found, i.e., q/qo ~ 7 - 10, the factor containing the function g 
and its derivative is of the order of the unity. Therefore, the en- 
ergy released by the reconnection is generally greater than the 
energy released by the helicity dissipation and the derivative of 
the surface term S in Eq. ( ITTb cannot be neglected. 

The chromospheric emission corresponding to the energy re- 
lease occurring in the corona is localized at the footprints of 
the flux tube that at each given time experiences a reconnection 
or a decrease of its helicity (cf. Sect. 12.6b . Therefore, a distant 
observer will see a modulation of the chromospheric emission 
with the orbital period of the planet since the amplitude of the 
modulation is independent of the orbital or rotation phase, given 
the axisymmetry of the field, and the visibility of the footprints 
is modulated by the motion of the planet along its orbit. This 
model prediction corresponds to the observations of the chro- 
mospheric hot spots i n HP 179949 and v And as reported by 
IShkolnik et all (I2005L 120081) . The phase lag between the longi- 
tude of the planet and the chromospheric footprints of the flux 
tube interacti ng with the planetary field has been discussed by 
lLanzal (120081) and we refer to that study for more details. 

2.9.2. Non-axisymmetric linear force-free field 

When |Z?i| is comparable or larger than |Z?ol, the situation is 
more involved. The power dissipated in the interaction when 
|Z?i| « |Z?ol is of the same order of magnitude of that released 
in the case of an axisymmetric field, but it depends explicitly on 
the time because (A p • B), 2B z A pz - B x A px , and R m are modu- 
lated by the terms sin 2 (0 - £) and cos 2 (0 - (cf. Eqs. [7] [33] 
and [34b. In other words, P lec and fhei are modulated with the 
frequency 2(a> - Q), i.e., twice the synodic frequency. From a 
rigorous point of view, the magnetospheric radius R m is a func- 
tion of B 2 which introduces contributions from the higher har- 
monics, i.e., 4(a) - Q), 6(a> - £2), etc., into the modulation of 
the energy dissipation rates. For simplicity, we assume that the 
fundamental frequency, i.e., 2(w - Q), dominates the variation 
of R m , leading to the same frequency for the modulation of P rec 
and Phei- Moreover, one should consider that the visibility of the 
chromospheric footprints in the case of a non-axisymmetric field 
is modulated by the rotation of the star which adds a further char- 
acteristic frequency, i.e. Q, and its harmonics to our variations. 
Typically, the first two harmonics, i.e., 20 and 3f2 should be 
considered. 

The combination of the modulation of P lec and Phei with the 
frequency 2(a> - Q.) and the modulation of the visibility of the 
footprints with the frequency Q and its harmonics produces a 



time dependence of the signal with frequencies (2a> - £2), (2a> - 
3£2), 2(l>, 2lo - 40, 2co - 50, and 2u> + O for a distant observer. 
Thus, the modulation of the chromospheric emission with the 
orbital frequency is no longer the dominant periodicity. This may 
lead to a confused situation in which it is virtually impossible to 
disentangle the SPMI signal from the intrinsic variability and 
the rotational modulation of the chrom ospheric emission typica l 
of an active star (see also the study bv lCranmer & Saa71l2007l) . 
This can explain why in several cases there has been no evidence 
of SPMI, even in those systems that have shown a modulation 
of the chromospheric emission with the orbital period in some 
seasons, such as HD 179949 or v And. Considering the Sun as a 
template, we see that the large scale magnetic field is close to be 
axisymmetric during the minima of the eleven-year cycle, while 
it shows a deviation from axisymmetry during the other phases 
of the cycle when large active regions or complexes of activity 
are present on the surface. 

The prediction that the modulation of the SPMI signal be- 
comes practically indistinguishable from intrinsic stellar activ- 
ity variations when the photospheric field has a sizable non- 
axisymmetric component can be tested by applying Zeeman 
Dopp ler Imaging techniq ues along the lines of the studies by, 



ppler imaging techniques along me lines ot i 
jMoutou et atl d2007l) and lFares et ail (1201 Ol) 



e-g 



2. 1 0. Non-linear fields 

In view of the high mathematical complexity of non-linear force- 
free fields, we limit ourselves to some specific kinds of axisym- 
metric fields that can be expressed as: 



B = 



rsinO 



r oO or R 



(37) 



where Bo sets the radial field intensity at the North pole, A(r, 6) 
is the flux function, and Q = Q{A) a scalar function that has a 
different functional form according to the specifically considered 
family of fields. Note that both A and Q are non-dim ensional in 
our de finition. The families of fields considered bv lLow & Lou I 
(1990) have a separable flux function of the form: 



A(r,0) = (r/R)- n f(x), 



(38) 



where x = cos 6, n is a positive constant, not necessarily an inte- 
ger, and / is given by the differential equation: 



1 



(1 - x l )f"(x) + n(n + l)f(x) + A'\1 + -j [/(x)] ,+//K - 0, (39) 

that is solved in [—1,1] subject to the boundary conditions 
/(— 1) = /(l) = with A 2 as an eigenvalue. The scalar func- 
tion Q in this case is given by: 



Q(A) = AA l+u ", 

while the force-free parameter a is: 

IdQ An+ l(ry l i 
a ~ ~R~dA ~ R n \r) 



(40) 



(41) 



We shall consider fields with < n < 1 because they have the 
slowest decay with the distance (cf. Eq. [381 and look for so- 
lutions of the boundary value problem for f(x) that satisfy the 
conditions f(-x) = f(x) in [—1,1] and /'(0) = 0, following the 
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method described bv lWolfsonl(ll995h . The magnetic field com 
ponents are: 

B r = -B a (rlR)- {n+2) f(x\ 

-(n+2) /(■*) 



B 6 = nB (r/Ry 



B = ABo(r/R) 



-(n+2) 



sin 8 
[/(*)] 
sinfl 



(42) 



The potential magnetic field with the same normal component 
2?p r = -B()f' ( x) on t he surface of the star is discussed in Sect. 5 
of IWolfsonl d 1995b . We consider only the component of its 
vector potential with the slowest decay with the distance from 
the star, i.e., that corresponding to the dipole component be- 
cause it dominates over the higher order multipoles given that 
the multipole of order k decays as (r/R)~ k . Therefore, we find 
A p = (3 / 4)BoR(r / R)~ 2 sin6>0 and assuming the planet in the 
equatorial plane: 

{A p -B) rp = ^B 2 RQ(A)^y . (43) 

To compute the energy dissipation rate from Eq. (13TV we con- 
sider the magnetospheric radius as given by Eq. ([7} and assume 
that the flux tube connecting the stellar surface with the bound- 
ary of the planetary magnetosphere has a nearly uniform force- 
free parameter, i.e., (a) ^ a, thus obtaining: 



Phel = £!l^ iaR f B ^/^ iA 2 +n 2 r u3 Vtel ( M- (n+5)/3 m 

16/zn + l p° p \RJ 

Since a is constant along magnetic field lines and is given by 
Eq. (|4TT >. we recast this equation as: 

_ 27 n n + 1 ,2^,^2,3^,2 , „2 r i/3 v / r Y ( " +U)l3 u ^ 

These equations are valid provided that a(A v ■ B) = aA p B^, is 
uniform over the surface of the planetary magnetosphere (cf. 
Sect. 12.6b . Since the relative variation of aBj, in the meridional 
direction increases with decreasing «, Eq. d45b is valid only when 
n ^ 0.1. For smaller values of n the energy dissipation rate is re- 
duced by a factor as nr/2R m , where r is the distance of the planet 
and R m the radius of its magnetosphere, as can be derived from 
Eq. ( f29b . Therefore, the energy dissipation rate does not diverge 
for n — > 0, but tends to a finite limit because the eigenvalue 
A 2 — > 1 in the same limit. The maximum value of the energy 
dissipation is obtained from Eq. (|43T > for n w 0.1. 

The magnetic field configuration s corresponding to our 
model have been studied in detail bv IWolfson I (11995 ). When 
the value of n decreases, the projections of the field lines on the 
meridional plane form taller and taller loops. In the limit n — > 0, 
the field lines become radially directed and an infinitely thin cur- 
rent sheet is formed on the equatorial plane. Such a configuration 
cannot support a confined corona because the hot plasma is free 
to escape along the radial field lines. For the existence of a con- 
fined corona, as indicated by the X-ray observations of planet- 
hosting stars, we should limit to values of n > 0.1. In this case, 
the total energy of the field is well below the so-called Aly limit, 
i.e., the minimum energy for the opening of all t he field lines for 
the assumed boundary conditions (see Sect. 4 in I Wolf son ll 19951 
for more details), and the pressure of the plasma cannot produce 
a loss of confinement of the corona by opening up its field lines. 

The ratio of the power released by the helicity modulation to 
that dissipated by the magnetic reconnection is: 



fhel 

p 

* rer 



27 1 n + 1 



16y r£ 



A 2 



A 2 I r V*" 1 
Tn 1 \R) ' 



(46) 



For a typical relative distance of the close-in planets (r/R) ~ 
7-10 and n ~ 0.1 the energy released by the helicity modula- 
tion is as 4 - 5 times larger than that produced by the reconnec- 
tion between the magnetic fields of the star and the planet. This 
implies that Eq. (l2"TT i is approximately valid for the considered 
magnetic fields and can be used to estimate the energy dissipa- 
tion rate neglecting the contribution of the time derivative of the 
surface term X. 

The footpoints of the field line connecting the stellar surface 
with the planet are located at colatitudes 9q and n-80, symmetric 
with respect to the equator. The value of 80 follows from the 
constancy of a along each field line and is the solution of the 
equation: 



/(COS ft>) = (!■/*)-", 



(47) 



since on the equatorial plane where the planet is located /(0) = 
1 . The azimuthal ang le Ac/> between th e footpoints and the planet 
is given by Eq. (8) of lWolfsonl (1 19951) that we reproduce here: 



n Jo 



[fix)] 
(1-x 2 ) 



dx. 



(48) 



Another class of axisymmetric force-fr ee fields of the kind 
specified by Eq. (l37l l has been proposed bv lFlver et al. I d2004). 
They have: 



[Q(A)f = ^jA k+ \ 



(49) 



where y is a parameter (constant for a given field configuration) 
and k an odd positive integer which ensures that the r.h.s. of this 
equation is positive definite independent of the sign of the flux 
function A. The force-free equation is solved subject to dipole- 
like boundary conditions, i.e., A(R, 8) - sin 2 8, |VA] — > as r — > 
00, and A(r, 8 = and n) = 0. With those boundary conditions, 
solutions to the force-free equatio n are found numerically only 
for k > 5 (see iFlver et ai~l l2004). Fields with k > 9 are very 
difficult to treat numerically owing to the steep gradients of Q. 

An interesting property of these fields is the possibility of 
an approximate asymptotic representation in the limi t r/R » 1 
by m eans of the Low & Lou's fields (cf. § 2.5 of IFlver et aTl 
2004). For k — 5, the representation by a field with n = 1/2 
is remarkably accurate, while for k — 7 and k — 9, the corre- 
sponding fields with n = 1/3 and n = 1/4, respectively, give less 
good representations in the range r/R ~ 7 - 10 of interest for our 
model. 



3. Results 

3.1. Coronal parameters and timescales 

We shall consider the systems for which some evidence of SPMI 
has been reported, notably HD 179949, v And, t Boo, and 
HD 189733. The stellar and planetary parameters adopted in our 
computations are reported in Table [T]where we list from the left 
to the right, the name of the system, the surface magnetic field 
intensity, the rotation period of the star, the orbital period, the 
semimajor axis of the orbit, the ratio a/R, and the relative veloc- 
ity v re i between the planet and the stellar coronal field. 

The magnetic fields of the host stars can be measured us- 
ing Zeeman Doppler Imaging techniques. For the F-type stars 
HD 179949, v And, and r Boo, we assume Bq = 10 G which is 
the m aximum radial field observed in t Boo close to the visibl e 
pole dCatala et alJl2007t iDonati et alJl2008t iFares et alJ 12009). 
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For the K-type star HP 1 89733 , we assume Bp = 40 G , as mea- 
sured by Mou tou et all (120071) and iFares et al.1 (l2010t) . the star 
being remarkably more active than the other targets. The rota- 
tion periods of HP 179949 v And, and HP 189733 are taken 
from IShkolnik et al l d2008l) . iPopp enhaeg er et al I (1201 ll) , and 
iHenrv & Winn! (120081) . respectively. ICatala et alfdioOTb found 
that t Boo has a remarkable differential rotation (as 20 percent) 
with a mean rotation period practically synchronized with the 
orbital motion of its planet. We assume a rotation period of 3.9 
days that corresponds to the high-latitude rotation period to max- 
imize the relative velocity between the coronal field and the or- 
bital motion of the planet. 

The X-ray coronal emissions of six planet hosting stars 
within 30 pc of the Sun have been fitted by Pop penhaeger et al71 
(2010) with a two-temperature model giving a cool component 
with a temperature of about 1 MK that dominates the emission 
measure. Only two stars show a comparable contribut ion from a 
hotter component with T ~ 4-5 MK. For HP 1 79949 iSaar et alj 
(2008) find a coronal temperature ranging from ~ 0.45 and 
~ 1.1 MK with the lower temperature component having an 
emission measure w 3 times la rger than that of the higher tem- 
perature component. In v And.lPoppenhae ger et al. I d201 ll) find 
a corona with a mean temperature of ~ 3 MK. 

We apply the coronal model of Sect. 12.11 to estimate the 
plasma B at the distance of the planets and the timescales char- 
acteristic of our systems in order to justify our approach. The 
variation of the magnetic field strength with the radial distance 
is comparable with that of a dipole field or is slower for our field 
configurations (Lanza 2009), thus we assume B(r) = Bp(r/R) . 
Indicating with Bp the value of the parameter at the base of the 
corona, we find: 

where Eq. (|4|i has been applied. Assuming an electron density 
«e0 = 10 14 irT 3 and a magnetic field of 10 G at the base of the 
corona, we find Bp = 0.007, 0.014, 0.027 for T = 1,2,4 MK, 
respectively. Considering a star with the same mass and ra- 
dius of the Sun, and a planetary distance r/R - 10, we find 
that B(r) = 10 -4 , 0.22, 13.8 for T = 1,2,4 MK, respectively. 
Therefore, we see that in a magnetostatic model the high tem- 
perature component of a stellar corona (T ^ 2 - 3 MK) cannot 
extend up to the distance of the planet because there B(r) > 1 . 
In other words, the hot component of the corona is confined into 
loops with a height of a few stellar radii, while the loops extend- 
ing up to the distance of the planet have a plasma temperature 
not exceeding 1-2 MK. 

The other crucial assumption of our model is that the 
timescale for the helicity variation tu is substantially longer than 
the relaxation time fR, that in turn is a few times the Alfven tran- 
sit time ?a along the magnetic field lines. The timescale fe is 
the time that the orbiting planet takes to travel across a distance 
equal to the diameter of its magnetosphere, i.e., fn = 2/? m /v re i, 
where v re i is listed in Table[T]for our systems. Assuming a typi- 
cal magnetospheric radius > 5 x 10 s m, we find ?h of the order 
of 10 4 s or longer for our systems. On the other hand, neglecting 
the curvature of the field lines, the Alfven transit time is approx- 
imately given by: 



Table 1. Stellar and planetary parameters 



(50) 



f 



dr 
VA(r) ' 



(51) 



Name 


So 


P rat 


Z'orb 


a 


a/R 


Vrel 




(G) 


(d) 


(d) 


(AU) 




(km s" 1 ) 


HD 179949 


10 


7.0 


3.09 


0.045 


7.72 


88.4 


HD 189733 


40 


11.9 


2.22 


0.031 


8.56 


125.4 


v And 


10 


9.5 


4.62 


0.059 


10.16 


71.4 


t Boo 


10 


3.9 


3.31 


0.046 


7.38 


22.8 



where v A (r) = B(r)/ -^up(r) is the Alfven velocity at the dis- 
tance r. Assuming that B(r) = Bp(r/R)~ 3 and adopting the model 



in Sect. 12.11 to compute the density p(r) with a base electron 
density n e0 = 10 14 nr 3 , B = 10 G, T = 1 MK, and a/R as 
listed in Tabled we find t A ~ 10 3 s for HP 179949 and r Boo; 
f A ~ 2.2 x 10 3 s for v And; and t A ~ 400 s for HP 189733 
because Bp — 40 G in that case. Note that the orbital veloc- 
ity of HP 189733 exceeds the isothermal sound speed for a 
coronal temperature < 10 6 K that can produce a hydrodynamic 
b ow shock in f r ont of the pl anetary magnetosphere according 
to lVidotto etal.1 (l2010bl 120111) . The shock is weak because the 
Mach number is ~ 1.1 - 1.2 so the magnetic field compression 
is negl igible for our purposes (see, e.g., § 5.2.3 of IVidotto et al.l 
2010a, for the effect of a strong perpendicular shock on the mag- 
netic field). Moreover, in the assumed low-beta environment of 
the stellar corona, the formation of a perpendicular shock re- 
quires that the relative orbital velocity of the planet exceeds the 
fast Alfven speed, i.e., V2/T 1/,2 c s , where c s is the isothermal 
sound speed, that is not verified in our case. 

In conclusion, the requirements on the timescales character- 
izing the processes considered in our model are generally satis- 
fied, provided that the temperature along the coronal field lines 
that touch the planetary magnetosphere is — 1 MK or lower, as 
required by the previous considerations on the confinement of 
the stellar corona. 



3.2. Linear force-free fields 

We compute the energy dissipation rates for the systems consid- 
ered in Sect. 13.11 Considering first an axisymmetric linear force- 
free magnetic field, we estimate the energy dissipation rates P KC 
and Phei from Eq. dSJ and Eq. d35l l. respectively, and compare 
the phase lag between the p lanet and the c hromospheric hot spot 
computed by the method of Lanza (2008) with the observations. 
Since we are interested in finding the maximum released power, 
we increase the value of a as much as possible because fhei 
scales as q^^q" 1 ^ that is |a| 4 ^ 3 since go = \cc\R and q - \a\a. 
However, a limitation on a is imposed by the requirement that 
the semimajor axis of the planetary orbit a be smaller than the 
limit radius = gi/M where g(q) has its first zero and the 
field lines close back onto the star. This implies in our cases 
\a\ < 0.5/T 1 . 

To compute the radius of the magnetosphere, for the non- 
transiting planets of HP 179949, v And, and r Boo, we assume 
a radius R p - 8.6 x 10 7 m, i.e., 1.2 Jupiter radii, while for 
the transiting planet of HP 189733 we have a measured radius 
R p = 8.08 x 10 7 m. The intensity of the planetary magnetic field 
at the poles is not known. Here, we adopt a field as strong as pos- 
sible, i.e., Bp = 100 G, which is about seven times the Jupiter's 

2/3 

field, because f rec and fhei scale as Z? ' and we aim at maxi- 
mizing the energy dissipation rate. Such extreme fields cannot 
be excluded in the framewo rk o f the dynamo models proposed 
bv IChristensen et ail (120091 ) and iReiners & Christensen (2010) 
provided that the internal heat flux of the planet, which controls 
the strength of the field, is sufficiently high. Since the dissipated 
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Fig. 2. Meridional section of the linear force-free magnetic field 
corresponding to the parameters assumed for HD 179949 in 
Table [2] Note the rope of azimuthal flux symmetric with respect 
to the equator and the position of the planet indicated by the 
filled dot on the equatorial plane (z - 0). The field line connect- 
ing the planet to the stellar surface is outside the flux rope while 
the field lines of the rope close onto themselves without reaching 
the stellar surface. 

2/3 

power scales as B L , adopting a field of 15 G, close to the value 
of Jupiter, will reduce the estimated powers by a factor of ~ 3.7. 

We choose the parameter bo/co defining the function g(q) 
to maximize the energy dissipation rate. The parameters of our 
field models are reported in Table [2]together with the energy dis- 
sipation rates. We list from the left to the right, the name of the 
system, the intensity of the stellar field Bq (see Sect. 13.11 for a 
justification of the chosen values), the intensity of the assumed 
planetary field B p q, the force-free parameter a, the parameter 
bo/ co, the footpoint colatitude 6q of the field line connecting the 
planet with the surface of the star, the phase lag A(p between 
the footpoints of this line and the planet, the dissipated power 
fhei as given by Eq. OSb . and the dissipated power P rec as given 
by Eq. ([8]). A meridional section of the magnetic field in the 
case of HD 179949 is plotted in Fig. [2] The field has a promi- 
nent azimuthal flux rope symmetric with respect to the equato- 
rial plane whose field lines are detached fro m the s t ellar s urface. 
This kind of configuration is discussed in lLanzal J2009). Note 
that the planet is located outside the flux rope, thus it is mag- 
netically connected with the stellar surface and can induce the 
formation of a chromospheric hot spot by releasing energy in 
the connecting flux tube. 



The dissipated powers listed in Table[2]can be compared with 
the observations of HD 179949 and u A nd that give powers o f 
* 10 2() and ss (2-3)x 10 19 W. respectively dShkolnik et al]200l . 
i.e., about two orders of magnitude larger than Phei and one order 
of magnitude larger than P lec as predicted by our model. Since 
the dissipated power scales as B^ 3 B^, to get P rec of the order of 

10 20 W, we should have Bq ~ 30-50 G, which is indeed the case 
for HD 189733 which is therefore the system with the highest 
predicted effect. The reason why there have been no conclusive 
observations of SPMI for this system may be its complex mag- 
netic field topolog y with a predominance of non-a xisymmetric 
components (cf. iMoutou et aTll2007t iFares et alJl2010h which 
lead to a complicated time dependence of the dissipated power 
(cf. Sect. |2.9.2| i. However, the possibility of a flaring activity 
modulated with the orbit of the planet in HD 189733 may indi- 
cate a remarkable SPMI in this system. 

Note that Phei < ^rec in all of our linear force-free mod- 
els. This indicates that the neglect of the surface term dLjdt in 
Eq. ( TP7I ). rigorously speaking, is not justified. The contribution 
of this term, however, is of the same order of magnitude of the 
power dissipated by the reconnection of the planetary and stellar 
fields, thus its inclusion is not expected to change our conclu- 
sions. 

An important constraint on the magnetic field model is pro- 
vided by the lag between the planet and t he phase of maximum 
visibility of the chromospheric hot spot. Shkol nik et al.l (120051) 
found that the maximum chromospheric emission falls at phase 
~ 0.7 - 0.8 in the case of HD 179949 and at phase- 0.5 in the 
case of v And, where phase corresponds to the inferior con- 
junction of the planet. In the case of t Boo, Walk er et al. 1 (l2008h 
found a pha se of maximum activity around ~ 0.7 - 0.8, while for 
HD 189733 IShkolnik et al.l (|2008) suggest a possible enhance- 
ment of the intra-night chromospheric variability with a phase 
lag of ss 0.8. These values correspond to angles of ~ 80° - 120° 
for HD 179949, r Boo, and HD 189733, and of ~ 180° in the 
case of v And. The A(p values derived from our models with 
\a\ = 0.5R~ l are at the upper bounds of those ranges for the first 
three stars, but are incompatible with the lag observed in v And. 
For this star, we cannot increase Iffl beyond 0.2/? -1 , otherwise it 
becomes impossible to reproduce the observed phase lag. As a 
consequence of the smaller \a\, we have a remarkable reduction 
of the dissipated powers. 

3.3. Non-linear force-free fields 

Here we limit o urselves to the axisymmetric fields discussed by 
Wolfson l (fl99l . A meridional section of the magnetic field lines 
in the case of HD 179949 is plotted in Fig. [3] for n = 0.5. The 
eigenvalue corresponding to n — 0.5 is A 2 = 0.82343. Note that 
this case gives also a good approxima tion to the more complex 
non-linear field of iFlver et al. I (120041) with k = 5 in the radial 
range of interest for star-planet interaction. In Table |3] we list, 
from the left to the right, the name of the system, the colatitude 
60 of the footpoint of the field line joining the stellar surface with 
the planet as given by Eq. d47l >. the azimuthal angle A<p between 
the footpoint and the planet as given by Eq. 0481 1. and the energy 
dissipation rates i\ei and P rec as derived from Eq. (|43T > and (|8), 
respectively, with Bq = 10 G in all the cases with the exception 
of HD 1 89733 for which B = 40 G. The planetary field strength 
BpO = 100 G in all the cases to maximize the energy dissipation 
rate as in the case of linear force-free fields. Assuming Z? p o = 
15 G - a much more realistic value in view of the models of 
iReiners & Christensen I (1201 Oh for planets with ages of at least 
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Table 2. Linear axisymmetric force-free fields 



Name 


B 


SpO 


\a\ 


bo/co 


8o 




-Phel 






(G) 


(G) 


(R- 1 ) 




(deg) 


(deg) 


(10 20 W) 


(10 20 W) 


HD 179949 


10 


100 


0.5 


-0.3 


40.0 


129.23 


0.09 


0.61 


HD 189733 


40 


100 


0.5 


-0.4 


20.0 


110.99 


0.68 


3.10 


v And 


10 


100 


0.2 


-1.1 


32.5 


177.1 


0.05 


0.01 


t Boo 


10 


100 


0.5 


-0.5 


32.0 


120.0 


0.05 


0.10 




Fig. 3. Meridional secti on of the non-linear force-free magnetic 
field of lWolfsonl d 19951) with n = 0.5 for the case of HD 179949. 
The filled dot indicates the close-in planet assumed to be on the 
equatorial plane of the star (z = 0) that is marked by the dashed 
line. 

Table 3. Non-linear axisymmetric force-free field for n = 0.5 



Name 


o 


Acp 


-Phel 


p 

1 rec 




(deg) 


(deg) 


(10 19 W) 


(10 19 W) 


HD 179949 


40.13 


53.99 


1.23 


0.44 


HD 189733 


38.66 


54.51 


1.04 


0.39 


v And 


36.28 


55.27 


0.35 


0.14 


t Boo 


40.78 


53.75 


0.38 


0.13 



1-2 Gyrs, and the field observed in Jupiter - the given powers 
are reduced by a factor of ~ 3.7. The variation of aB^ across 
the planetary magnetosphere does not exceed 5 percent in all the 
cases, thus our assumption of a constant coronal field over the 
volume of the planetary magnetosphere is well justified. 

For n = 0.25 {A 2 = 1.01203), P h d is a factor of 2.4 larger, 
but is still insufficient to account for the observations, while the 



ratio Phei/Aec ranges from 2.8 to 3.5. However, n = 0.25 gives a 
better agreement with the observed phase lags, because in all the 
cases A(f> ranges between 68? 1 and 69°.7 with 6q between 57°9 
and 6Q°.l. We can decrease n up to, say, n = 0.1, because for this 
value the variation of aB^, across the planetary magnetosphere 
reaches ~ 60 percent making our formula only roughly valid. 
However, even in this extreme case, Phel increases by a factor of 
4.4 with respect to the case with n = 0.5 and this is in the best 

case only marginally compatible with the observations; 

We conclude that the non-linear force-free field o f lWolfsonl 
(119951) cannot account for the dissipated powers observed in our 
systems, even for extreme values of the planetary magnetic field. 
For HD 179949, r Boo, and HD 189733, it can account for the 
phase lag between the planet and the chromospheric hot spot if 
n = 0.25, but in the case of v And the predicted lag is too small. 
An insufficient dissipated power is expe cted also in t he ca se of 
the more general non-linear models of Fi ver et al. I (2004) be- 
cause the Low & Lou models give a fairly good approximation 
to them in the radial range considered for the star-planet interac- 
tion. 

4. Discussion and conclusions 

We have introduced an approximate model to compute the en- 
ergy dissipated in the interaction between the magnetic field of 
a star and a close-in planet. Our model assumes that the orbital 
velocity of the planet is much smaller than the Alfven speed in 
the stellar corona and that the coronal field perturbed by the mo- 
tion of the planet relaxes to its minimum energy state within 
a timescale comparable with the Alfven transit time along the 
coronal field lines. With these assumptions, we can treat the evo- 
lution of the field as a sequence of magnetostatic configurations 
and estimate the energy variations under the constraint that the 
total magnetic helicity of the field is conserved. 

We have estimated the power dissipated by magnetic recon- 
nection as well as that released by the modulation of the field he- 
licity associated with the orbital motion of the planet. The latter 
process can operate also when the planetary field is very small or 
even absent because what is needed is a relative velocity between 
the planet and the coronal field of the star. If the planetary field 
BpO = 0, we can assume that the radius of the planetary magneto- 
sphere corresponds to the radius of the planet, i.e., R m = R p and 
apply, e.g., Eq. (03]l simply substituting B^^B^R 2 with B 2 R 2 . 

When a planetary dipolar field is present, we find that the 
dissipated power scales as R 2 B^ 3 B 2 J Q 3 , with Bo the surface field 
of the star. This is a general result, independent of the specific 
mechanism responsible for the magnetic energy dissipation, pro- 
vided that the dissipated power is proportional to the available 
magnetic energy, that scales as B 2 r and the surface of the plane- 
tary magnetosphere, that scales as B Q 2/3 B 2 J Q 3 R 2 . This scaling law 
may be used to infer the relative planetary field strength, not yet 
directly observable, from the measurement of the power of the 
SPMI in a sample of stars, as suggested bv lScharf I (120101) . 
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An important conclusion of our model is that, while the re- 
connection of the stellar and planetary fields releases an addi- 
tional power in the corona that can lead to an excess X-ray emis- 
sion, this is not the case for the energy released by the helic- 
ity variation that produces a modulation of the dissipation of 
the energy already available for the coronal heating or flaring. 
Since the power released by the reconnection is generally of 
th e order of 10 17 - 10 1 8 W, t his can expla in the recent results 
of Canto Martins et aT71 (1201 ll) or those of [Poppenhaeger et al. I 
that suggest at the most a marginal correlation of the stel- 
lar X-ray emission with the mass and the inverse orbital semi- 
major axis of the planet. On the other hand, a modulation of the 
co ronal flaring activity w ith the orbit of the planet, as suggested 
bylMllitteri et aT| d201 lb, is predicted by o ur model. Concerning 
the correlation found by lHartman I (|20 1 Ob . the range spanned by 
the stellar chromospheric emission versus the planetary surface 
gravity covers approximately one order of magnitude. This ap- 
pears to be too large to be accounted for by the energy released 
by the interactions considered in our model and requires a dif- 
ferent explanation. 

The time modulation of the chromospheric or X-ray emis- 
sion predicted by our model has a single frequency, equal to 
the orbital frequency of the planet, only if the stellar magnetic 
field is predominantly axisymmetric, as expected when the star 
is close to the minimum of its activity cycle by analogy with the 
Sun. When the stellar field has a non-axisymmetric component 
comparable or larger than the axisymmetric one, the modulation 
becomes multiperiodic with frequencies coming from the com- 
bination of the orbital and stellar rotation frequencies and their 
harmonics, which makes the modulation virtually indistinguish- 
able from the intrinsic activity fluctuations of the star. This may 
explain why in some seasons the chromospheric signature of the 
SPMI has not been observed in stars that had previously shown 
some evidence of modulation with the orbital period. 

Models of the SPMI based on stellar linear force-free fields 
allow us to treat both the case of axisymmetric and non- 
axisymmetric fields with an analytical description. However, the 
power released, even assuming the extreme values allowed for 
the free parameters, is insufficient by at least one order of mag- 
nitude. This is related to the existence of an upper bound for the 
force-free parameter a because the field must extend at least up 
to the radius of the planetary orbit. On the other hand, such con- 
figurations can be useful to describe the final state of a confined 
magnetic structure, when its excess magnetic energy has been 
dissipated and the field has reached the minimum allowed en- 
ergy compatible with the conservation of the relative magnetic 
helicity. Such closed field configurations have been considered 
by, e.g jLanza d2010t) when discussing the evolution of the rota- 
tion of stars with hot Jupiters. 

Non-linear force-free fields in general do not have limita- 
tions on the value of the force-free parameter a. However, the 
classes of fields that we have considered provide dissipated pow- 
ers that are still insufficient by at least one order of magnitude to 
account for the observations, although they can account for the 
phase lag between the planet and the hot spot, with the exception 
of v And. Therefore, our results cast doubts on a straightforward 
association of the observed chromospheric hot spots with the or- 
biting planets, requiring an alternative mechanism t o explain the 
phenom enon. Indeed, the numerical simulations of Cohenetal] 
(1201 lbl) . including also the kinetic energy of the plasma flows in 
the system, are more promising and could be more appropriate 
than the present idealized models to explain the observations. 
In those numerical models, there are magnetic loops intercon- 
necting the star with the planet. Therefore, the maximum power 



made available by the relative motion of the planet, as derived 
from the flux of the Poynting vector [T l E x B across the base 
nRp of an interconnecting flux tube, where E — -v x B is the 
electric field and v the relative velocity, is: P nfX^R^B^v. 

For a relative velocity v = 10 4 - 10 5 ms~',a planetary field 
BpO = 10 G, and a radius of the planet R v — 7 x 10 7 m, we obtain 
a maximum available power of ~ 10 - 10 21 W that could be 
enoug h to account for the hot spots observed by Sh kolnik et alJ 
(2005). We shall explore this possibility in more detail in a forth- 
coming work. 
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Appendix A: The vector potential of a potential 
magnetic field 

The application of Eq. )3 U requires the determination of the vector potential A v 
for a given distribution of the normal component of the magnetic field B n on the 
surface of the star. We consider a spherical polar coordinate system (r, 6, <f>) with 
the origin at the barycentre of the star and the polar axis along the rotation axis 
of the star. The surface of the star is a sphere of radius R. 

The normal component of the magnetic field at the surface is B,-(R, 8, <f>) and 
is equal to the radial component of the potential magnetic field B p = Vt/r, i.e., 
(di///dr)n = B r (R, 9, <j>). The general expression of the scalar potential of the field 



B„ exterior to the star is: 



<Kr,8,<f>) = ^^[a,,, 
n=l m=0 



,YL t (6,<f>) + b nm r„ 



„<ft#>](J) 



-C*+l) 



(A.l) 



where a nm and b nm are the numerical coefficients of the expansion of the 
scalar potential in terms of the spherical harmonic functions Y^ m (8, cp) = 
P™(cos8)cos(m<j>) and 1^(9,0) = P^(cos 9) sm(m0), and P% are Legendre as- 
sociated polynomials of degree n and azimuthal order m with < m < n. The 
term corresponding to n = is absent in the development because the magnetic 
field has no monopoles. If the field is interior to the sphere of radius R, Eq. 1A. U 
is still valid provided that the radial dependence of the component of the order n 
is changed as (r/R)". 

Thanks to the orthogonality properties of the spherical harmonics, we can 
immediately derive the coefficients a,„„ and b„,„ from a surface integration of the 
radial field component, i.e.: 



j f*1ll r*K 

Nmn JO JO 



RB r (R,e,4>)Y% 



)sm8d8d<p, 



Y r-ln r>n 
"inn Jo JO 



8, <f>)Y;; m (0, <j>) sinO d8d<f>, 



where the normalization factor 



- An 



(n +1) (n + m)\ 
(2n + 1) (n-m)!' 



(A.2) 



(A.3) 



The connection between the scalar potential ip(r, 8, <j>) and the vector potent ial 
A p can be derived as follows. In the representation of Chandrasekharl J 19611) . a 
potential field has only a poloidal component, i.e., it can be written as: B p = V X 
jV X [<t>(r, 8, 0)r]|, where r is the unit vector in the radial direction and <D(r, 6, <j>) a 
scalar function that will be specified below. Therefore, the corresponding vector 
potential satisfying V ■ A„ = is: 



V X [<D(r, 8, <f>)r] ■■ 



1 



rsin6\d(f>J 



—\e- -(—^' 



06 j T ' 



(A.4) 



where 9 and are the unit vectors in the meridional and the azimuthal directions, 
respectively. To derive the func tion <P, we compare the components of the field 
Bp as given by the formulae of Chandrasekhar (1961) with the components of 
Vifi. Chandrasekhar's formulae are: 



B pr = i/Jo, B p 



1 d ld®\ 



rdeXdr)' Bf *~' 



i a /<9<d 



r&md d<f> \ dr 



where 
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L~ = — — | sine- 
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Since 



is the angular part of the Laplacian operator, i.e., V 2 = ^ { r2 §?) 
the field B p is potential, V 2 ip = 0. From this equation and the first of Eq. jA.51 . 
it is easy to see that if/ = 3? . The same expression immediately verifies also the 
second and the third relationships in JA.51 ; hence we find: 



<D(r, 8. 



<p) = J if/(r',8,</>)dr' = 

°° 1 " Y ~ 



(A.7) 



where a„„, and b nm are determined by the normal component of the magnetic 
field at the surface of the star according to Eqs. dA.21 . By substituing Eq. )A.7) 
into Eq. IA.4K we derive the components of the vector potential A p . Note that 
when the field is axisymmetric (d/d</> = 0), A p has only the azimuthal compo- 
nent, i. e., it is everywhere parallel to the orbital velocity of the planet assumed 
to move on an equatorial and circular orbit. 



Appendix B: Perturbation of the energy of the 
coronal field 

The presence of the planetary magnetosphere produces a perturbation of the en- 
ergy of the coronal field. If the planet were absent, the total energy of the field 
would be: 



Eo 



X 



B z dV, 



(B.l) 



where B is the coronal field and V the volume exterior to the star. When the 
planet is present, the coronal field is B m = B + B', where B' = and the po- 
tential f is given by Eq. )24t . Therefore, the energy of the coronal field becomes: 



1 f 



2H 



Bt dV, 



(B.2) 



where V is the volume of the corona exterior to the star and the planetary mag- 



netosphere. Since Bj. 



B + 2B ■ B' + B = B + (B + B') ■ B' + B ■ B' 



B 1 + (B + B') ■ Vf + B ■ V£ = B 1 + B m ■ Vf + B ■ Vf = B 1 + V ■ (fB m ) + V ■ ((B), 
we can apply Gauss' theorem to find: 



2fi{L 



B z dV + 



<r(B m ■ h) ds + 



£(B ■ h) dS\ 



(B.3) 



where S a is the surface of the star and 5 ra the surface of the magnetosphere. The 
potential f = (1/2)(B ■ rXRj^/r 2 ) is negligible on the surface of the star, while 
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on the boundary of the magnetosphere B m ■ h = 0. Therefore, the middle integral 
vanishes and in the third integral only the integration over S m remains, yielding: 

El = ±- f B 2 dV - -!-« m £ (B ■ ff dS, (B.4) 

2/J Jy, Ajl J Sm 

where the minus sign in front of the second integral comes from the orientation 
of the normal to S m , i.e., h = -r. Adopting a reference frame with the z axis 
along the relative velocity v between the planet and the coronal field, and the 
xz plane containing the vector B, as specified in the text, we can perform the 
integration and find: 

El = — I B 2 dV - — V m B 2 , (B.5) 
2/J Jv 4[i 

where V m = (4ff/3)i?^ is the volume of the magnetosphere. Since the radius of 
the magnetosphere is given by Eq. 0, we recast this equation as: 



~ f 



B dV - — RpBpoB, (B.6) 

where R p is the radius of the planet and B p o the planetary field at the pole. Finally, 
the available energy difference between the coronal field configurations with and 
without the planet is: 



3 t n 



AE = Ei - E = - — V m B' = —SiBjoB, (B.7) 
An fx 

where we made use of V = V U V m , with V m the volume of the magnetosphere, 
and the energy of the coronal field in that volume when the planet is absent is 
jy^(B 2 /2fi) dV = V m B 2 /2/i, provided that B can be assumed uniform over the 
volume V* ra . 



Appendix C: Linear force-free field and vector 
potential of the corresponding potential field 

The components of the lin e ar for ce-free field of order n = 1 as introduced by 
Chandrasekhar & Kendall] d!957h can be written as (see Appendix [A] for the 
adopted reference frame): 

Br = (®) g(9)[B o cos6> + fi 1 sinecos(0-f)l, 

1 <?n 

Bo = --— |g'(?)[B O sine-B 1 cos0cos(^-f)] + 

2 q 

+ g(a)BisH<t>-0} 

^ 2 

B4, = -— [ g (q)B sm8 + Big'(q)sin(<p-O+ 
2 q 

- fli£te)cos0cos(0-£)], (C.l) 

where Bo is the intensity of the axisymmetric mode of the field, i.e., that with 
m = 0, at the poles and Bi is the intensity of the radial component of the non- 
axisymmetric mode with m = 1 at the longitude ^. Note that the axisymmetric 
field has a zero radial component on the equatorial plane because the field is 
similar to that of a dipole; the dimensionless radial distances q = \a\r and qo = 
\a\R, while the function g is defined as: 

g( q) = qZl(ci) = [ fc o J -3/2(g) + c J3, 2 (q)] V? (C 2) 

qoZliqo) [br,J-3/2(qo) + Co/3/2(?0)] V?0 ' 

where Z| is defined by Eq. (10) of Chandrasekhar & Kendall 1 1957) for n = 
1, J-3/2 and J3/2 are Bessel functions of the first kind of order -3/2 and 3/2, 
respectively, and b p and cp are coefficients determined by the field at the base 
of the corona (see Lanza 2009, for more details); the function g' = dg/dq is 
the first derivative of g. When comparing the pres ent expression of the magnetic 
field with that of Chandrasekhar & Kendall 1 1957), note that in the second line 
of their Eq. ( 1 3) there is a typo, thus the correct argument of their radial derivative 
is [rZ„(ar)]. 

The components of the vector potential of the potential magnetic field hav- 
ing the same radial component of the force-free field at the surface, computed 
with the method introduced in Appendix lAl are: 

Ap r = 0, 

A pS = -i«Bisin(<S-#f^ 
-V = jR[Bosine-Bicosf?cos(^-f)]|^J . (C.3) 
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